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Abstract 

We discuss the binding energy Ej, of impurities in semiconductors within density functional 
theory (DFT) and the GW approximation, focusing on donors in nanowires as an example. We 
show that DFT succeeds in the calculation of Ef, from the Kohn-Sham (KS) hamiltonian of the 
ionized impurity, but fails in the calculation of Ef, from the KS hamiltonian of the neutral impurity, 
as it misses most of the interaction of the bound electron with the surface polarization charges of 
the donor. We trace this deficiency back to the lack of screened exchange in the present functionals. 
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The binding energy Eh of donors and acceptors is a key quantity in semiconductor physics 
because it determines the doping efficiency In semiconductor nanostructures for example, 
confinement and electrostatics tend to shift the impurity levels deeper in the gap, which 
decreases dopant activity jl, 2]. Therefore, ab initio calculations of impurity binding energies 
are highly desirable to assess the performances of ultimate nanodevices. Besides, donors 
and acceptors are the prototypes of charged defects in semiconductors, and a fundamental 
understanding of the strengths and weaknesses of present ab initio approaches such as density 
functional theory (DFT) and the GW approximation {3], 4] would open the way to a more 
accurate modeling of complex defects. 

So far, the calculation of Eh in bulk semiconductors has been possible only with semi- 

nn 

empirical methods [5|, [6[ . However, calculations based on DFT have become practicable in 
ultimate nanostructures with a smaller number of atoms. Recently, the case of donors in Si 
nanowires (Si NWs) has been adressed with both semi-empirical methods and DFT, with 
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contradictory results. Tight-binding |7|, |8] and effective mass calculations [9!], supported 
by experiments [l, 2], indeed suggest that Eh increases as 1/R with decreasing wire radius 
R, due to the interaction of the bound electron with the surface polarization (or "image") 
charges of the impurity, resulting in a significant decrease of the doping efficiency in the 



R < 10 nm range. In contrast, DFT calculations lOj, lllf predict that Eh decreases much 
faster than 1/R, and is about 3 — 4 times lower than found in Refs. [7, 9]. In this Letter, 
we show that present DFT approaches, based on the Kohn-Sham (KS) hamiltonian of the 
neutral donor, can not predict Eh correctly in bulk and nanostructures, because they miss 
most of the interactions of the carriers with the polarization charges of the impurity. We 
propose an alternative strategy based on the KS hamiltonian of the ionized donor which 
circumvents this deficiency. 

For a donor, Ef, is the energy needed to ionize the neutral impurity and bring the electron 
to the conduction band edge far away. It can be defined as the difference Eh = I d (N + 1) — 
A P (N) between the ionization energy I d (N+l) of the neutral impurity (with N+l electrons) 
and the affinity A P (N) of the pristine system (with N electrons and no dopant). Since 
A(N) = I(N + 1), the binding energy can also be computed as an isoelectronic difference of 
ionization energies, E b = I d (N+l)-I p (N+l), or affinities, E b = A d (N)-A p (N). In practice, 
the ionization energies and affinities can be computed either as total energy differences 
I(N) = E(N-1)-E(N) and A(N) = E(N) - E(N + 1) [10], or as "quasiparticle" energies 
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Uj , i.e., as the highest occupied (HOMO) and lowest unoccupied (LUMO) molecular orbital 



energies. However, the quasiparticle problem should in principle be adressed with many- 



body perturbation theories (MBP' 



to miss the HOMO-LUMO gap 



s) such as the GW approximation, since DFT is known 



13| . We actually show hereafter that the above 



definitions of Eb are consistent in the GW approximation, but not in DFT. Using the insight 
gained from many-body theory, we conjecture that DFT should succeed in the calculation 
of Eb from the KS LUMO of the ionized impurity, but fails in the calculation of Eb from 
the KS HOMO of the neutral impurity, due to the lack of explicit screened exchange in the 
present functionals. We support these conclusions with DFT calculations on Si NWs. 

The binding energy in many-body theory - In MBPT, the quasiparticle energies E n and 
wave functions (p n of the iV-electron system are the solutions of the quasiparticle equation: 

-^A r <^ n (r) + v ion (r)if n (r) + v h (r)ip n (r) 

+ J dV S xc (r, r', E n )y n {r') = E n cp n (r) (1) 

where fi on (r) is the ionic potential, t>h(r) = / d 3 r'p(r')/|r — r'| is the Hartree potential 
created by the ground-state electronic density p(r), and S xc (r, r', E n ) is the "self-energy" 
that describes exchange and correlation effects. The ionization energy is I(N) = —En, while 
the affinity is A(N) = —E N+1 . 

The GW approximation has become the reference for the calculation of the band structure 



of semiconductors 
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For illustrative purposes, we shall use hereafter the simpler static 



COHSEX form (COulomb Hole and Screened EXchange) of the GW self-energy p]. 
then be split in two parts Ecoh + ^sex'- 



£ CO H(r,r') = ^5W N (r, r)5(r - r') (2a) 
SsEx(r,r') = -WMr,r') 7 (r,r') (2b) 
7 (r,r') = V*Ar) Vn ,(r'), (2c) 

n'e{(r} 

where the sum runs over the occupied states with a given spin a. H^r, r') is the screened 
Coulomb interaction, i.e. the total potential created at point r' by a test unit charge at 
point r 14]. It can also be split in the bare potential v(r, r') = l/|r — r'| created by this test 
charge, plus the response 5Wn(t, r') = Wn(t,t') — v(r,r') of the valence electrons. Esex 
has the same functional form as the Hartree-Fock exchange, but with a screened instead of a 



bare Coulomb interaction. EcoH( r ; r describes the interaction of a carrier at point r with 
the valence electrons which dynamically respond to its motion. 

Before adressing the impurity problem, we shall discuss the form of Wn(t,t') in bulk 
materials and nanowires. In a solid, a test charge q t — +1 at point r attracts valence 
electrons in a small "cloud" around (over ~ a bond length). This cloud contains a total 
charge q c = —(1 — 1/k), where k is the static dielectric constant of the material. The 
electrons are actually dragged from the surface of the system, where they leave an opposite 
polarization (or "image") charge q s = —q c . In bulk, these image charges are infinitely far 
away, so that the long-range potential created by the test charge is simply W^(r,r') ~ 
(q t + <7 C ) /|r — r'| = 1/ (/t|r — r'|). In a nanowire, however, the electrons are dragged within a 
few -R's only from q t , so that the transfer of charges from the surface to the cloud becomes 
shorter and shorter-ranged with decreasing R. The screening is therefore reduced by q s and 
the potential ultimately tends to Wn(t, r') ~ l/|r — r'| when R — > (i.e., the test charge 
mostly sees vacuum around for small -R's). This simple picture is consistent with classical 
electrostatics (where the surface polarization charges are given bv the discontinuity of the 
electric field), and fully supported by quantum calculations la, Il6 |. 

As discussed previously, the many-body binding energy of an impurity can be computed 
as Eb = A d (N) — A P (N), the difference between the affinities of the ionized impurity and 
pristine systems. They fulfill the equation H p ' d (N) V p / +1 = —A p ' d (N){p p T f +1 , where H P (N) 
and H d (N) are the respective quasiparticle hamiltonians: 

H p ' d (N) = -~ A + + vf + Sg x + EgJ H . (3) 

The physics of the impurity is most easily brought out from the difference between H P (N) 
and H d (N). On one hand, the extra proton of the ionized impurity is screened by the valence 
electrons through the Hartree potential v£- Neglecting short range chemical corrections in 



a first approximation 13] , we can therefore write: 

Kn + <]-Kn + <]^-^(r„r), (4) 

where is the impurity position. On the other hand, we do not expect significant differences 
between the screened Coulomb interactions Wfj and W^, nor between the one-particle den- 
sity matrices 7 d and 7 P , except possibly right around the donor and surface (image charges), 
on length scales much shorter than the Bohr radius of the impurity. Hence, Sg EX ~ Sg EX , 



yd ^ y'P nnrl- 
^COH — ^COH' ana - 

H d (N) ~H P (N)- W N (r h r). (5) 

In a first approximation, the quasiparticle hamiltonian of the ionized impurity is the quasi- 
particle hamiltonian of the pristine system plus the screened Coulomb potential of a unit 
charge at the impurity position. This is the usual "hydrogenic model" Q| used in Refs. jl-9| 
to calculate Eb in Si NWs. 

The electron is therefore bound to the impurity by the screened Coulomb interaction 
Wjv( r i) r )- I 11 bulk silicon, Wjsr^i^v) ~ l/[/«|r^ — r|] and Eb ~ 50 meV. In a nanowire, 
however, the electron also interacts with the image charges of the donor. Since the total 
surface polarization charge is q s = (1 — 1/ k) 3> 1/k, this leads to a large oc 1/R enhancement 
of Eb with decreasing R p]]. 

Let us now compute the binding energy E b = I d (N + 1) - P(N + 1) from the ionization 
energy of the neutral impurity. I (N + 1) and I P (N + 1) fulfill the equation H p,d (N + 
1)^jv+i = —I p ' d (N+ l){p P jf +1 , where, as before, H P (N+1) and H d (N+l) are the quasiparticle 
hamiltonians of the (A r +l)-electron pristine and impurity systems. In the latter, the HOMO 
<Pn+i is the occupied bound state of the impurity. The neutral impurity as a whole now 
introduces a localized perturbation of the pristine system which is screened by the valence 
electrons. We can therefore write: 



v. 



where: 



+ <] - K n + <\ * -WZ+iin, r) + vt(r) , (6) 

<r) = J dV< +1 (r,r')|^ +1 (r')| 2 (7) 

accounts for the screening of the bound state potential. Assuming again that W^ +1 ~ W^ +1 , 
and that the valence band wave functions (pi,.--,<PN are little affected by the neutral impurity, 
we further get: 

S^ EX (r, r') - S p SEX (r, r') ~ -W N+1 (r, r') 

x[< +1 (r)^ +1 (r') - < +1 (r)^ +1 (r')] . (8) 

The second term can be neglected in bulk and nanowires where <f p N+1 is an extended state. 
The first term cancels f&(r) when applied to the HOMO ^jv+i- The effective hamiltonian 
for the bound electron therefore reads: 

H d (N+l)~H p (N + l)-W N+1 (r h r). (9) 



In principle, I(N + 1) = A(N) and we should have recovered the same equation as before 
[Eq. flS}]. Here H P (N) is however replaced with H P (N + 1) and W N with W N+ i. Since 
^iv+i i s an extended state, H P (N + 1) and H P (N) also primarily differ by the substitution 
Wn — > Wn+i- The appearance of W^+i introduces a residual "self-correlation" error in 



the GW ionization energies [18j], which is however expected to be limited in solids. We 
can therefore conclude that GW provides a consistent description of the binding energies, 
whether computed from A d (N) or I d (N + 1). 

This paragraph clearly demonstrates the importance of screened exchange in the calcu- 
lation of I d (N + 1). Screened exchange indeed cancels the unphysical screened interaction 
of the bound electron with itself which arises from u&(r) in Eq. H d (N + 1) is there- 

fore the hamiltonian of a charged system as expected (the bound electron interacts with 
N + 1 ionic charges but N electrons). Such spurious self- interactions are a serious issue in 
self-consistent descriptions of occupied localized states. In this respect, we would like to 
point out that the Hartree-Fock (HF) bare exchange £ x (r, r') = — v(r, r')7(r, r') does not 
properly correct the screened self-interactions appearing in solids. Following the same lines 
as before, the HF hamiltonian of the HOMO of the neutral impurity can indeed be written 
H*p(N + 1) ~ H^ F (N + 1) - W N+l (r u r) + <(r), where: 

<(r) = J d 3 r' [WWr.rO -«(r,rO] l^+i(r')| 2 - (10) 

f^ r (r) is the spurious potential created by the valence electrons in response to the bound 
state density |<y9^ +1 (r)| 2 , i.e. the potential created by a diffuse charge p e ff( r ) — (1 — 1/k) x 
| V^a^+i ( r ) | 2 pl us the opposite surface polarization charge q s = —(1 — 1/k). These surface 
polarization charges balance those embedded in the impurity potential Wjv+i( r i) r ) (equiv- 
alent, as discussed before, to the potential of a net charge 1/k at rj and q s — (1 — 1/k) at 
the surface). H^ F (N + 1) is therefore approximately equal to the hamiltonian of the pristine 
system plus the bare Coulomb potential of a unit charge spread around the impurity (the 
charge 1/k at the impurity position plus the diffuse charge p e ff around) . As a consequence, 
Peff plays the role in the HF approximation of an effective polarization charge, mislocalized 
within the scale of the Bohr radius instead of the surface. The relative error on Ej, should be 
limited in thin nanowires (where the Bohr radius is comparable to R) , and maximum in bulk. 
The implications for hybrid functionals in DFT will be discussed in the next paragraph. 
The binding energy in DFT - We now discuss the binding energy of the donor within 
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DFT. For the sake of simplicity, we first focus on the local density (LDA) and generalized 
gradients (GGA) approximations, then address the case of hybrid functionals. In DFT, the 
ground-state density p(r) of the iV-electron system is computed from the eigenstates of the 



Kohn-Sham hamiltonian 



19|: 



- ^A r Lp n (r) + [v ion + v h + w xc ](r)y? n (r) = E n Lp n (r) (11) 

where v xc (r) is the exchange-correlation potential. In LDA and GGA, v xc (r) = v xc (p(r)) is 
a function of the local density p(r) and of its derivatives. DFT is known to underestimate 



the band gap energy of semiconductors [13j. Still, we show below that DFT should succeed 
in the calculation of the binding energy from the LUMO of the ionized impurity, but that 
present functionals fail on the neutral impurity. 

Let us first compute = A d (N) — A P (N) from the LUMOs of the KS hamiltonians 
H^ S (N) and H^_ S (N). The previous arguments are also valid in DFT: The extra proton 
of the donor is screened by the valence electrons, so that Eq. (j4j) still holds. We do not, 
moreover, expect much differences between v p c (r) and v xc (r), except possibly right around 
the impurity. Therefore, in a first approximation: 

HiaW-HSsW-WNfar). (12) 

We hence recover the hydrogenic model as before [Eq. (jSJ)]. The KS hamiltonian of the 
ionized impurity thus embeds the same extra physics (with respect to the hamiltonian of 
the pristine system) as the GW approximation: Although the LUMO energies are typically 
underestimated by DFT, the binding energies computed as the difference between the KS 
LUMOs of the ionized impurity and pristine systems should be reasonably accurate. This 
only holds, of course, as long as the binding energy is not too large with respect to the DFT 
band gap. 

Let us now compute Ej, = I d (N + 1) — I P (N+ 1) from the HOMOs of the KS hamiltonians 
H^ S (N + 1) and H^ S (N + 1). The KS wave function <p% +1 is the occupied bound state of 
the impurity. The neutral impurity as a whole is again screened by the valence electrons 
[Eq. (j§J)]. The exchange-correlation potential v xc (r) is also affected by the extra bound state 
density around the impurity. We hence get : 

H^{N + 1) ~ H^ s (N + l)-W N+1 (r h r) 

+ v b (r) + Av xc (r) , (13) 



where At> xc (r) = v d c (r) — v xc( r )- At variance with the GW approximation, At> xc (r), a local 
density correction within the Bohr radius, can not be expected to cancel t>b(r) [Eq. ([7j)], a 
long-range Coulomb term. This results in i) a self-interaction error, and ii), an almost com- 
plete cancellation of the interaction of the electron with the image charges of the impurity. 
Indeed, Wjv+i( r i, r) and Vb(r) are the potentials created by two opposite charges (the ionized 
impurity and bound electron), leaving no net charge in the hamiltonian. Both errors give 
rise to an increase of the impurity level and to a decrease of the binding energy. Although 
this is especially sensitive in thin nanowires, where the enhancement of Ef, is mostly due to 
the interaction with the image charges, the LDA and GGA would fail up to the bulk [where 
the impurity potential decreases exponentially instead of l/(ft|r — r^)]. We stress that the 
calculation of the ionization energy or affinity of the impurity as a difference of total energies 
[lo| . which involves the neutral impurity as the initial or final state, suffers from the same 
deficiencies in the LDA or GGA. 

Application to Si nanowires - The binding energies of dopant impurities in Si NWs have 



been previously computec 
hybrid functional (HGH) 



from the KS HOMO of the neutral impurity using GGA and a 



111 ], i.e. a mixture of Hartree-Fock bare exchange with GGA. As 



discussed previously, bare echange does not localize the polarization charges properly, the 
error being however likely limited in thin nanowires (the total charge being correct). The 
GGA results of Ref. [ll] are therefore expected to completely miss image charge effects, while 
the HGH results, which include 12% bare exchange, are expected to account for ~ 12% of 
the interactions with the image charges (even though mislocalized). As a consequence, the 
difference between the GGA and HGH results of Ref. should be approximately 12% of 
the image charge correction given by Eq. (3) of Ref. that is, 0.12 eV for R = 1 nm, 0.17 
eV for R = 0.75 nm, and 0.25 eV for R = 0.5 nm. This is actually in good agreement with 
the data of Table I of Ref. [11 ]. 

To further support the above conclusions, we have computed the binding energy of a P 
impurity at the center of a hydrogen passivated, [110]-oriented Si nanowire with diameter 
d = 1.73 nm, either as Ej, = I d (N + 1) - A P (N), or as E£ = A d (N) - A P (N), using KS 
HOMOs and LUMOs as ionization energies and affinities. The LDA was used in a wavelet 
basis set as implemented in the BigDFT code [20]. The neutral impurity was first relaxed in 
a 660 atoms supercell. Since the treatment of a charged system is still problematic within 
such a supercell approach, A d (N) (as well as I d (N + 1) and A P (N) for consistency) were 
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actually computed in finite rods with lengths I up to 9.2 nm (1584 atoms). These rods were 
build from the original 660 atoms supercell by connecting segments of pristine nanowires 
and hydrogen passivated ends. The binding energies computed from the charged and neutral 
impurities are respectively = 0.93 eV and E[ = 0.06 eV for I = 9.2 nm. As expected, 
E^ is much larger than E[, and in good agreement with the semi-empirical model of Ref. [t| 
(Eb = 0.92 eV when / — > oo). This confirms that present functionals are able to predict the 
binding energies of impurities or defects from the KS hamiltonian of the charged defect. 

To conclude, we have shown, by a formal comparison with the GW approximation, that 
the donor binding energies computed from the Kohn-Sham hamiltonians of neutral impu- 
rities can be strongly underestimated in semiconductor nanostructures (even with hybrid 
functionals). This is due to the lack of screened exchange in the present functionals, and 
explains the discrepancies between Refs. Q, and previous works 0, 0]. The binding 
energy of a donor should preferably be computed as the difference between the Kohn-Sham 
LUMOs of the ionized impurity and pristine systems. This provides a reasonable substitute 
for much more expensive GW calculations of defect bound states in solids. 
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